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. . , Abstract 

OO 

0> ' In the framework of the perturbative QCD and the diquark model of baryons we 

^^ . have obtained the fragmentation functions for heavy quark that split into polarized 

Afe baryons. We predict the longitudinal polarization asymmetry for the prompt A^ 

produced in e~^e~ annihilation at the Z resonance and estimate the spin-1/2 and 

spin-3/2 beauty baryon production rate. 






O 

f^ '. Introduction 

> . 

^ ' The heavy baryon production in e+e^ annihilation is an increasingly important subject 

^r^ ■ to study. The recent measurements near the threshold of the heavy quark production, as 

C^ . well as at the Z resonance, give us the information about heavy baryon masses, life times, 

decay modes Jl], ^, || and polarizations @, |^. 

0\ [ The b quarks produced in e~^e~ annihilation at the Z peak are strongly polarized. 

^ ' Accordingly to the Standard Model the right-left longitudinal asymmetry is A^^j^ = —0.94. 

'P' This polarization is expected to be only slightly (2%) reduced by the hard gluon 

Q-i. emission before the hadronization into baryons. Therefore the measurement of the A;, 

r-| ! polarization is very important for study the details of the hadronization processes. The 

first measurement ^ gave an intriguingly small value Aj^j^ = — 0.23lo'23) which contradicts 



X 



the prediction that the large part of the initial polarization of the quark to transfer to the 
Afo baryon. 
^ The data on the A;, polarization may be described in part (A^*^ ^ —0.68) by means of 

a simple depolarization model, based on the cascade A;, partial production in the decays 
of higher spin beauty baryon produced in the Z decays [0]. 

In this paper we propose the depolarizing mechanism occurring during the prompt 
Aj, production via b quark fragmentation. Our approach base on the perturbative QCD, 
the quark-diquark model of the baryons [Q and the nonrelativistic approximation, which 
is used successfully for the description of the large distance effects in the heavy quarko- 
nium production H. The same approach have been discussed in the our previous paper 



10], where the predictions for the doubly heavy baryon production via the heavy quark 



fragmentation have been presented. 



1 The prompt A^ production in Z decays 

The fragmentation function Dh^j>^^ {z, fi) at the initial scale fi = fio = frtq + 2m£) for the 
production of the A^ baryon, containing the Q = b quark and the scalar diquark D = (ud), 



is given by the following expression ||T0|, pJ]] : 
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where Ai is the matrix element for the production Aj, baryon and antidiquark D with the 
total four-momentum q = p+q' and invariant mass s = g^ in the Z decay {Z —>■ AhDb), Aio 
is the matrix element for the production of a real b quark with the same three-momentum 
q. The lower limit in the integral (1) is 
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where M = itlq + m£, is the baryon mass, mg is the heavy quark mass, tud is the diquark 
mass. 

The amplitude A^o may be presented as follows: 



Mo{Z ^ hi) = U{q)T''V{q)e^{Z), 
where ea{Z) is the polarization four- vector of the Z boson and 

r" = - _ ^^„ 7" (ct. - C^7^ 



(2) 
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is the quark-boson vertex. 

In the axial gauge for the gluon propagator associated with the four-vector n 
(1,0,0,-1): 



daf3{k) = -gaf3 + 



(kn) 



the fragmentation contribution comes only from the Feynman diagram shown in Fig.l. 
After some obvious simplification we have obtained the matrix element of the Z decay 
into Afe baryon ||10|| : 
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2U{p) ' 
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2 \2 



r"\/(g)£,(Z) 



(3) 



where g'^ = 4:Tcas, 45'-' /3v^ is the color factor of the diagram, Fs{k'^) is the form factor 
of the vertex g* -^ DD, ^^(O) is the A;, = (QD) wave function at the origin in the quark- 
diquark approximation. The gluon coupling to scalar diquarks was used in the following 
form 0: 

S^^ = -ig,T'^{q'-pn).F,ik'), (4) 

where T" = A"/2 are Gell-Mann matrices, k = q' + p^, is the gluon four-momentum, the 
spinor Lf{p) describes A;, baryon, p = Pg + Pn, Pd = rp, pg = {1 — r)p, r = mo/M. The 
scalar diquark form factor at fc^ > may be parameterized as in ref . [|12| , where the authors 



fit successfully the angular distributions of the baryons in the processes 77 
the widths of J/%Ij — > pp, AA: 



pp, AA and 



(5) 



where Q^ = 3.22 GeV^ and form factor is restricted to value smaller than 1.3. We use 
also the parameterization (5) with the fixed full width at half maximum V '^ 0.8 GeV^. 
In the both cases there are no singularities in the physical region of the virtuality of the 
fragmenting quark or the square of the gluon four-momentum. 

The phenomenological diquark form factor parameterization (5) corresponds to the 
elastic vertex g* -^ DD [|l^. However, in the fragmentation processes the contribution 
from the inelastic vertex g* -^ Dud may be dominant. There are no any information 
about inelastic diquark form factors. That is why we will use parameterization (5) to 
describe spin effects in the A?, production. As it will be shown, the choice of the form 
factor parameterization is important only for the prediction the absolute values of the 
heavy baryon production rates, but the spin asymmetry doesn't depend on a diquark 
form factor in the kinematic region, which is studied. 

The reduced mass of the heavy quark - diquark system is the same order as one for the 
system of the two charmed quarks and we can hope that the calculation of the parameter 
\l/(0) using the potential approach, will be well-grounded and may be done with the quark- 
diquark interaction potential, which has been fixed previously in calculating the heavy 
baryon mass spectrum |13| . 

The energy distribution of the right (R) or the left (L) longitudinally polarized A;, 
baryons produced in the Z decays reduces at leading order in a^ to 



where 






P0 + P3 



(6) 



qo + qs 

p = {po, 0, 0,^3) is 4-momentum of the Af, , q = (go, 0, 0, gs) is 4-momentum of the frag- 
menting b-quark and /i ~ 'mz/2. Let us define 



The unpolarized fragmentation function is 

L>6^A,(2;,/i) = DI^_^j^^{z, fi) + D^_^^^{z, fi). 
The longitudinal asymmetry Ap^i may be presented as follows: 



A^l(^,/i) 



In the limit of r7i| >> s, M^ we have got 



AD 



b^hi 



{z,fx) 



D^ 



b^Ah 



>,/^) 



(7) 



(8) 



(9) 



Tr 



gr„gT^] {-g-^ + Z^Z^/ml) = mU&^ + C^^^) 



9^ 
4 cos^ 9w 



(10) 



Tr 



■g75r«tr J i-g'^f + Z'^Z^I 



Ulr, 



2 (~\b (~\b 



-2miCtC 



g' 



4 cos2 e 



w 



'11^ 



Taking into account that the four-vector of the A^ polarization in the states with the 
different longitudinal polarization is Sf^^ = ±(p3/M, 0, 0,po/^) and using the following 
exact formulae 



2pk = s — rriQ, k'^ = r{s — mg), 



2pg = s-mg + 2M^(1 -r), 2Zg = m| + s - m^, 

we can reduce the squared matrix elements to the form, which contains expressions (10) 
or (11). 

Omitting the next details, we write here the obtained results for the unpolarized 
fragmentation function at the scale n = fio'- 
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z^ 
(1-^(1 -r))2' ^^^^ 

The result for the ADb^\^ {z, Ho) may be written as for the Db^A^ {z, jio) but we don't 
present it here because of it's unwieldy. 

The recalculating of the fragmentation functions D^^{z, fio) from the start point fi = 
fio to the fi > fio may be done using Gribov-Lipatov-Altareli-Parisi (GLAP) evolution 
equations 

dD , , [^ dy _ ( z ^ 
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where Pq^q is the splitting function at leading order in ag' 
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f{x)+ = f{x) - 6{1 -x) f{x')dx'. 

Jo 

Note that at leading order in as one has 

J PQ^Qiz,fx)dz = 0, 
and the evolution equation implies that the fragmentation probability 

Pb^At = / Db^At{z,fJ^)dz 



as well as the total asymmetry Aj:j^ don't evolve with the scale /i. However, the 2;— dependence 
of the Aj^i^i^z, n) changes drastically when the scale /z growth from /i = /Zq to /i ~ razjl- 

First of all, we discuss here the results concerning spin effects in the prompt A^ produc- 
tion, which don't depend on the value of \E'(0). Figs. 2-4 demonstrate the dependence of 
the polarization asymmetry A^*^ on the diquark mass, the form factor parameterization, 
z and QCD evolution scale \i. 

We can see in Fig. 2 that A^^ 7^ A^^^ ^^ the realistic values of the diquark mass 
(m£) = 0.6 — 0.9 GeV corresponds to \P^^i\ = 0.90 — 0.87) and only at the limit of 
m£) -^ one has Aj^j^ ^ Aj^^ = —0.94. Our result corrects the assumption, which is 
used usually, that the A^ baryon retain a large part of the initial polarization of the b 
quark [jl4|. The predictions for the asymmetry as a function of z strongly depends on the 
diquark mass at /i = /Xq, but at the large /i this dependence vanishes (Fig. 3). 

Apl'i{z,^) may be sensitive to the choice of the form factor parameterization at the 
small /i (see Fig. 4), however at the scale /i = mz/2 the results corresponding the different 
form factor parameterizations are equal. The measurement of the A^*^ versus z at the 
different /i is a exact test of the hadronization model. However, it may be done exper- 
imentally if we can to separate the prompt A^ production in the b quark fragmentation 
from the A;, produced in the cascades: 6 ^ S^, ... plus the hadronical decays S^, ... — * Af,. 
This opportunity was shown recently in the study of the heavy quarkonium production 
at FNAL []T3|, where the vertex detector technique was used. 

The Figs. 5 show our results for the Ajf^ versus z at the different /i = /Xq, mz/2 and 
rriD = 0.6 GeV. We have obtained A^"^ = —0.58 opposite to A^^ = —0.67. The probabil- 
ity of the heavy quark spin turning during the hadronization is about 13% for the c-quark 
and only 4% for the b-quark. The both values are independent from the diquark form 
factor and were obtained at niD = 0.6 GeV. We see that the simple rule (~ mj^/mo) is 
satisfied. 

2 The cascade A^ production in the Z decays 

Using the value of the radial part of the Af, baryon wave function |i?(0)p = 0.8 GeV^ (a^ = 
0.4, 771]:, = 0.6 GeV), which was calculated in the quark-diquark approximation |T^, we 
have found the total probability of the b quark fragmentation into the Af,: Pb^At ~ 
2.5 • 10~^ in the case of the form factor (5) and Pb^A^ ~ 2.5 ■ 10~^ if the singularities 
are restricted using the fixed width F. In spite of this deference the shapes of the frag- 
mentation functions Dh^j^^lz, fi) are the same for the both form factors (Fig. 6). The 
obtained value Pb^A^ gives the number of A^ per the Z decay which is smaller than it was 
measured at LEP ^. We conclude that the main part of the Af, baryons are produced in 
the cascade processes via the decays of the higher spin b baryons. The DELPHI collabo- 
ration has presented evidence for S^ baryon production |1^ and measured the rate to be 



Pb-^T.* = 4.8 ± 1.6%, the total b baryon production rate to be Pb^baryon = 11.5 ± 4.0% 0. 
We can estimate the fragmentation probability for the b quark that split into S^ 
baryon in the our approach. However, as it was demonstrated in the case of the Af, baryon 
production, the predicted absolute production rate is strongly depends on diquark form 
factor. It will be more reliable to calculate the ratio of the fragmentation probabilities for 



the b quark that spht into vector diquark states to scalar diquark states, using deferent 
form factors. As it was shown in ref. this parameter (called A) controls the A(, 
polarization in the cascade production. 

The gluon coupling to axial vector diquarks is presented by the following expression: 



V^ = tgsT''^e},eDiq'-pn),F,ie) - [(g'es)^D^ - (Pd^d)^5m]^2(A:' 



ie*DQ'){eDPD){q'-PD)^.F3{e 
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where e*jy,e£) are the diquark polarization four- vectors. Fi,F2 and F3 are form factors 
depending on the momentum transfer squared k"^ = {q' + PoY- The same as in ref. [^ 
we put for /c^ > 0: 
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where Qv = 1.50 GeV^ and k being the anomalous chromomagnetic moment of the vector 
diquark. The anomalous magnetic moment of the (ud) vector diquark was determined 
in ref. ]T2|], but it is not obviously that the chromomagnetic and magnetic moments are 
equal. We will choose k = for the simplicity. 

The amplitude for heavy quark fragmentation into spin-3/2 baryons, corresponding to 
fusion of the heavy quark and the axial vector diquark, can be written as a sum of two 
parts, which are proportional to Fi and F2 form factors: 
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Here, the ^fj.{p) is the Rarita-Schwinger spinor for the spin-3/2 baryons. The procedure 
of the calculation for the fragmentation function Dh^Y,*{z, fi) is the same as for the A^ 
baryons. 



Using the set of parameters: |-R(0)| = O.SGeV 



as 



0.4, niD = 0.6 GeV, /t = 0, 



we have found that the ratio A of the fragmentation probabilities for the b quark that 
split into vector diquark states to scalar diquark states is approximately equal to 5 ± 1 



independently on the choice of the diquark form factor. Accordingly the formula (5.10) 
ref.[|7i one has 



A, l + {l + AwMI^ 



K'l = \^/ ' Kl, (22) 

where the parameter wi gives the probability that spin-1 diquark has maximum an- 
gular momentum j^ = dzl along the fragmentation axis. The authors ref. M used 



A = 0.45, Wi = 0, as it is motivated by the Lund fragmentation model [|I^, and ob- 
tained that A^*^ = —0.68. If we put A = 5 in (22), as it follows from our results, we find 
that A^*^ = —0.24 which is in a good agreement with the value of the measurement. 

In conclusion we want note that the prediction of the cascade model for the A^*^ 
in the Z decays was obtained using the assumption that the turning of the b quark spin 
may be ignored, as is appropriate to the heavy quark approximation. Our results show 
that the probability of the b quark overturning during the hadronization may equal to 
4-7 % and the more careful analysis of the depolarization in the cascade A;, production is 
needed. 

When our work was completed the result of the computing the heavy baryon pro- 
duction rates at CLEO and LEP energies, which was obtained a similar way, has been 
presented [^ . 
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Figure captions 

1. Diagram used for description of the process Z° -^ Kj)D. 

2. The asymmetry A^*^ as a function of the ratio m£,/M. 

3. The asymmetry Apl'^ as a function oi z at fi = fio (curves 1 and 2) and /i = mz/'2 
(curves 3 and 4). The solid lines correspond to mr, = 0.9 GeV, the dotted lines 
correspond to tti/j = 0.6 GeV. 

4. The asymmetry A^*^ as a function of z at rriD = 0.6 GeV and fi = fio (curves 1 
and 2), mz/2 (curves 3 and 4). The solid lines correspond to diquark form factor 
parameterization (5), the dotted lines correspond to formula (5) with the fixed width 
r = 0.8 GeV2. 

5. The asymmetry Ajf^ as a function of z at rriD = 0.6 and jjl = jjlo (curve 1), mz/2 
(curve 2). 

6. The fragmentation function Db^A^{z, fi) normalized to unity at /i = /^o (curves 1 
and 2) and /i = mz/2 (curves 3 and 4). The solid and dotted lines correspond to 
the different form factor parameterizations. 



Figure 1: Diagram used for description of the process Z° — * Ky})D. 



Figure 2: The asymmetry K^j^ as a function of the ratio itid/M . 
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Figure 3: The asymmetry A^*^ as a function of 2; at /i = /io (curves 1 and 2) and /i = 
m^/2 (curves 3 and 4). TThe solid lines correspond to niD = 0.9 GeV, the dotted lines 
correspond to itld = 0.6 GeV. 
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Figure 4: The asymmetry A^*^ as a function of z at mD = 0.6 GeV and n = fio (curves 
1 and 2), rnz/^ (curves 3 and 4). The sohd hues correspond to diquark form factor 
parameterization (5), the dotted hues correspond to formula (5) with the fixed width 
r = 0.8 GeV2. 



12 



Figure 5: The asymmetry A^j^ as a function of z at m/) = 0.6 and [i = [lo (curve 1), 
mz/'2, (curve 2). 
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Figure 6: The fragmentation function Dh-tAi,{z, /i) normalized to unity at /i = /xq (curves 
1 and 2) and /i = mz/l (curves 3 and 4). The sohd and dotted hues correspond to the 
different form factor parameterizations. 



14 



